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Abstract. This paper concerns with the convergence analysis of a fourth order singular per- 
turbation of the Dirichlet Monge-Ampere problem in the n-dimensional radial symmetric case. A 
detailed study of the fourth order problem is presented. In particular, various a priori estimates with 
explicit dependence on the perturbation parameter e are derived, and a crucial convexity property 
is also proved for the solution of the fourth order problem. Using these estimates and the convexity 
property, we prove that the solution of the perturbed problem converges uniformly and compactly to 
the unique convex viscosity solution of the Dirichlet Monge-Ampere problem. Rates of convergence 
in the i?*-norm for fc = 0, 1, 2 are established, and illustrating numerical experiment results are also 
presented in the paper. 
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1. Introduction. In this paper we consider the following fourth order singularly 
perturbed Monge-Ampere problem: 

-eA^M"^ + det{D^u^) = / in n, (1.1a) 

(1.1b) 
(1.1c) 

where e > 0, C R" is a convex domain and / > in fi. The above problem is a 
special case of the following fourth order singular perturbation problem: 

eA^u^ +F(D2^i,V^i^^i^x) = in r2, (1.2a) 

on dil, (1.2b) 
Au' = £ ondn (1.2c) 

with F{D'^u'',\Ju'',u'',x) = f{x) - det{D'^u''). Here, D'^u^{x) and Wu^{x) denote 
respectively the Hessian matrix and the gradient of at a; e il. 

Problem (1.1 1, respectively (1.2), was proposed by the authors in as a step- 
stone to develop efhcient numerical methods, particularly Galerkin-type methods, 
for approximating the viscosity solution of the Dirichlet Monge-Ampere problem (cf. 

mm)- 

det{D^u) = f in 17, (1.3a) 
u = g on dft, (1.3b) 

respectively, the general fully nonlinear second order PDE problem: 

F{D^u, Vm, u,x)^0 in n, (1.4a) 
u^g on dn. (1.4b) 
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We remark that since the boundary condition (1.1c I (resp., (1.2c)) is artificially im- 
posed to make the perturbed problem well-posed; other type boundary conditions can 
be used in the place of ( l.lc[ ) (resp., (1.2c)) (cf. [H]). 

Fully nonlinear second order PDEs, which have experienced extensive analytical 
developments in the past thirty years (cf. [IHIIB]), arise from many scientific and en- 
gineering fields including differential geometry, optimal control, mass transportation, 
geostrophic fiuid, meteorology, and general relativity (cf. [71 [HJ [TSl [21 [Ml [31 [H] and 
the references therein). Fully nonlinear PDEs play a critical role in the solutions of 
these application problems because they appear one way or another in the governing 
equations. By the middle eighties the classical solution theory for fully nonlinear sec- 
ond order PDEs was well established (cf. [16, Chapter 17]). The classical framework 
was soon followed by a comprehensive viscosity (weak) solution theory (cf. [21 [7] and 
the references therein) after the introduction of viscosity solutions by Crandall and 
Lions in 8J for fully nonlinear first order Hamilton- Jacobi equations. 

In contrast with the success of PDE analysis, numerical solutions for general fully 
nonlinear second order PDEs was a relatively untouched area until very recently (cf. 
[TT] and the references therein). The lack of progress is mainly due to the following 
two facts: (i) the notion of viscosity solutions is nonvariational; (ii) the conditional 
uniqueness (i.e., uniqueness only holds in a restrictive function class) of viscosity 
solutions is difficult to handle at the discrete level. The first difficulty prevents a 
direct construction of Galerkin-type methods and forces one to use indirect approaches 
as done in |171 1121 114] for approximating viscosity solutions; this is main focus of 
this paper. The second difficulty prevents any straightforward construction of finite 
difference methods because such a method does not have a mechanism to enforce the 
conditional uniqueness and often fails to capture the sought-after viscosity solution 
(cf. [TT] and the reference therein). To overcome the above difficulties, inspired by 
the vanishing viscosity method for the Hamilton- Jacobi equations, in |12j we proposed 
to approximate the fully nonlinear problem (1.3) (resp. (1.4)) by the fourth order 
quasilinear problem (1.1) (resp. (1.2)). We called this approximation procedure the 
vanishing moment method. If the vanishing moment method converges, then the 
solution of the former can be computed via the solution of the latter. This can be 
achieved by using many existing numerical methodologies including all Galerkin-type 
methods. So far, all of the numerical experiments and the analysis of [1211131 fl4l[23l [22] 
show that the vanishing moment method works effectively. 

The primary goal of this paper is to provide a rigorous proof of the convergence of 
the vanishing moment method applied to the Monge- Ampere equation in the radially 
symmetric case in dimension n > 2. Specifically, we shall give a detailed analysis of 
problem ( |1.1[ ) in the radially symmetric case and prove the convergence of the solution 
of problem ( |1.1[ ) to the unique convex viscosity solution u of problem (1.3). We 
also derive the rates of convergence for u — u"^ in various of norms. These results then 
put the vanishing moment method on a solid footing and provides a partial theoretical 
foundation for the numerical work in [12l [131 IHl ES] . 

The remainder of the paper is organized as follows. In Section [2] we give the 
radially symmetric descriptions of problems (1.3) and (1.1). Sections |3] and |4j which 
are the bulk of this paper, is devoted to studying the existence, uniqueness, regularity 
and convexity of the solution of the singularly perturbed problem ( [1.1[ ) . With these 
results in hand, we prove the uniform convergence of the solution of problem ( [1.1[ ) 
to the unique convex viscosity solution u of problem (1.3) in Section [sj In addition, 
we derive rates of convergence for u — in the _ff*''-norni for A: = 0, 1, 2 in Section [6| 
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We present a few illustrating numerical computational results in Section [7] and end 
the paper with some concluding remarks in Section [8j 

2. Descriptions of radially symmetric problems. Standard space notation 
is adopted in this paper, we refer the reader to [H [TB] for their exact definitions. 
Unless stated otherwise, fl = i3/j(0) C R" (n > 2) stands for the ball centered at 
the origin with radius R. We do not assume O is the unit ball because many of our 
results will depend on the size of the radius R. The unlabeled letter C is used to 
denote a generic positive constant independent of e that may take on different values 
at different occurrences. 

Suppose that / = /(?"),/ ^ and g = g{r) in (1.3), that is, / and g are radial. 
Then the solution u of (1.3) is expected to be radial, namely, u{x) is a function of 

r := \x\ — \J'Y^^=i ^'j- We set u{r) :— u{\x\) — u{x), and for the reader's convenience, 

we now compute Au, A^m and det(D^u) in terms of ii (cf. [lHHS]). First, a direct 

calculation shows — Xj/r and (^^^^^ = By the chain rule we have 

— = Ur(r)- — = Ur{r) — , - — — ^ - -Ur{r) XiXj H dij. 

oxj ox J r dxjOXi r\r Jr r 

Here, the subscripts stand for the derivatives with respect to the subscript variables. 

Noting that D'^u{x) is a diagonal perturbation of a scaled rank-one matrix xx , 
and since the eigenvalues of xx'^ are (with multiplicity n — 1) and |xp = (with 
multiplicity 1 and corresponding eigenvector x), the eigenvalues of D^u{x) are Ai :— 

"''^-''^ lur('') ) = Urr{r) (with multiplicity 1), and A2 := -^^^^ (with multiplicity 



1). Thus, 

n — 1 1 

Au{x) = Ai + (n - 1)A2 = iirrir) H —Ur{r) = -^—^{r^'^Ur)^, 



Ur{r) 



n-1 I 



det{D^u{x)) = Ai(A2)"-l = iirrir) [-^\ = {{UrDr- 

Abusing the notion by denoting u{r) by w(r), problem ( |1.3[ ) becomes seeking a 
function u = u{r) such that 

^(K)"V = / iii(0,i?), (2.1a) 
u{R)^g{R), (2.1b) 
Ur{Q) = 0. (2.1c) 



We remark that boundary condition (2.1c) is due to the symmetry of it = u{r). 

Lemma 2.1. Suppose that r"-'^f e L^{{0,R)) and f > a.e. on (0,i?). Then 
there exists exactly one real solution if n is odd and there are exactly two real solutions 
if n is even, to the boundary value problem ( |2.1[ ). Moreover, the solutions are given 
by the formula 

g{R) ± (nL f{s)) " ds if n is even, 

u{r) = { ' ^ (2.2) 

g{R) ~ J^ i^nL f (s)^ " ds if n is odd 
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for r e (0, i?). Here, the function Lf{s) is defined by 



Lf{s) := / f'-'f{t)dt. 
Jo 



(2.3) 



Since the proof is elementary (cf. [211 \'25\). we omit it. Clearly, when n is even, 
the first solution (with the "+" sign) is concave, and the second solution (with the 
"— " sign) is convex because Ur and simultaneously positive and negative 

respectively in the two cases. When n is odd, the real solution is convex. 

Remark 2.1. The above theorem shows that u is at a point ro G (0, R) as long 
as f is at ro and Lf(rQ) ^ 0. Also, u is smooth in (0,i?) if f is smooth in (0,i?). 
We refer the reader to [211 I25j for the precise conditions on f at r — Q to ensure 
the regularity of u at r = 0, extensions to the complex Monge- Ampere equation, and 
generalized Monge-Ampere equations in which f = f{\/u,u,x). 

Similarly, it is expected that = M^(r) is also radial, and the vanishing moment 



approximation (1.2) then becomes 



u'iR)^g{R), 

;(0) = 0, K^(0)|<cx), \u^„^{r)\ ^ o{-^) asr^ 



ul^{R) 



n-1 
R 



uliR) 



(2.4b) 
(2.4c) 

(2.4d) 



In the next sections, we shall analyze problem ( |2.4[ ) which includes proving the ex- 
istence, uniqueness and regularity of the solution. After this is done, we then show 



(2.4c) and limr_j.o+ r*^ 



that the solution u'^ of (2.4) converges to the unique convex solution of (2.1). 

Integrating over (0,r) after multiplying (2.4a) by r"~^, using boundary condition 



-7/'- — "^^ i;*^ 



er 



— 0, we get 



(2.5) 



Introduce the new function w^{r) := r" u^{r). A direct calculation shows that 
w'^ satisfies 



CT"-if^<) + ^(w^)" = Lf in (0,i?). 



Converting the boundary conditions (2.4c)-(2.4d) to w'^ we have 



w'^(O) =<(0) =0, 



In addition, since — r" ^u^,. + (rt — l)r" ^u^, and w:;.^ — r" "m: 



l)r"- ^w^^ + (n — l)(n — 2)r" ■'u^, we have 



(0) = of for < j < min{2, n - I). 



(2.6a) 

(2.6b) 
(2.6c) 

2{n - 
(2.7) 



In summary, we have derived from (2.4a) a reduced equation (2.6a), which is 



second order and easier to handle. After problem (2.6) is fully understood, we then 



come back to analyze problem (2.4) 



3. Existence, uniqueness, and regularity of vanishing moment approx- 



imations. We now prove that problem (2.6) possesses a unique nonnegative classical 



solution. First, we state and prove the following uniqueness result. 



Theorem 3.1. Problem (2.6) has at most one nonnegative classical solution. 



Proof. Suppose that wf and are two nonnegative classical solutions to (2.6a) 



(2.6c). Let 



and 



a+fS=n-l 

a,p>a 



— Wo 



otherwise. 



Subtracting the corresponding equations satisfied by wf and yields 



-er 



'-(0) 



--€{0)^0^ 



(3.1a) 

(3.1b) 
(3.1c) 



Noting that > in [0, i?], we have by the weak maximum principle |10( Theorem 
2, page 329] max[o.«] 10^01 = max{|0-(O)|, 10^^)1} = niax{0, W{R)\}. If (t>'{R) = 
0, then (jf = 0. If <jf{R) ^ 0, then 0^ takes its maximum or minimum value at 
r = R. However, the strong maximum principle |241 Theorem 4, page 7] implies that 
(priR) 7^ 0, which contradicts with boundary condition (t>%{R) = 0. Hence, <j)^ = 
which implies wf = D 



Next, we prove that the existence of nonnegative solutions to problem (2.6) 



Theorem 3.2. Suppose r'^'^f e L^{{0,R)) andf>0 a.e. m (0,i?). Then there 



is a nonnegative classical solution to problem (2.6) 



Proof. We divide the proof into three steps. 

Step 1: Let 1/'° € C^([0,i?]) be a nonnegative function that satisfies = 
'0r(O) = and ipri-^) = ei?"~^. For example, we may take tp^{r) — f?'". We then 
define a sequence of functions {ip''}k>o recursively by solving for fc = 0, 1, 2, • • • 



-er 



nr"^ 

V;'=+l(0) = V.^+1(0) 

i/,^+l(i?) =£i?' 



0, 



in (0,i?), 



(3.2a) 

(3.2b) 
(3.2c) 



We first show by induction that for any such sequence satisfying (3.2), there holds 



> in [0, R] for all fc > 0. Note that 1/)° > by construction. Suppose that ip'^ > 
in [0, R]. Since / > and / ^ 0, we have 



fe+1 



> 



in (0,i?). 



Hence, ip'^'^^ is a supersolution to the linear differential operator appearing in the 



left-hand side of (3.2a). By the weak maximum principle [101 Theorem 2, page 329], 
we have 

minV'''+^(r) > TJim{0,iP''+\0), iJ/+\R)} = mm{0, ij/+^R)}. 

[0,-R] 



Now suppose that Tp''+^{R) < 0. Then since £^""^ > 0, the strong 

maximum principle [24, Theorem 4, page 7] imphes that ■0*^+-'^ = const in [0,i?], 
which leads to a contradiction as ip''~^^{Q) = 0. Thus, we must have ip'^'^^{R) > 0, and 
therefore, ■0'^+^ > in [0,i?]. By the induction argument, we conclude that '0*' ^ 
in [0, R] for all fc > 0. 

It then follows from the standard theory for linear elliptic equations (cf. [I0| 116)) 
that (3.2) has a unique classical solution 0*^+^. Hence the {k + l)th iterate tp'^^^ is 
well defined and therefore so is the sequence {'0'^}fc>o- 

Step 2: Next, we shall derive some uniform (in fc) estimates for the sequence 
{0'^}/c>o- To this end, we first prove that 0'^"'"^(i?) can be bounded from above 
uniformly in fc. Multiplying (3.2a) by ■0'^+^ and integrating by parts yields 



r=R 



£0«+^(r)0'^-+^(r) +e / Wl ^^ 

r=0 



(3.3) 



n- 1 
r 



i)''+\r)Lf{r) dr. 







It follows from boundary conditions (3.2b) and (3.2c I that 



e0^'+i(r)0^+i(r) ' ^= -e^i?"-! 0fc+i(i?). 

r=0 



Integrating by parts gives us 



(3.4) 



-0,^+1 (r)0'=+i (r) dr = ^ (0^=+! (r)) ' 
r 2r ^ ' 



2r2 



(0'=+i(r))'dr 



1 

2R 



2r2 



{tP^+^{r)f dr. 



(3.5) 



By Schwarz, Poincare, and Young's inequalities, we get 



il)''+^{r)L}{r) dr < 



B. ^2 rR 

H'r^'irrdr+'g 



{Lf{r)fdr (3.6) 



for some positive constant Ci — Ci{R). Combining (3.3|-(3.6| we obtain 



2e2i?"-i0'=+i(i?)+e / \^^+'{r)\Ur + 



e(n- 1) 



R 



(0'=+i(r))^dr+ / _^(0'=(r))"-i|0'=+'(r)pdr 



< 



nr 
2 rR 



{Lf{r)ydr. 



(3.7) 



Let z 0''+i(i?), b := f^, and c := j^(^(L/(i?)) . Then from we have 
— 62: — c < 0, which in turn implies that 



b ~ + 4c 6 + + 4c 
■Zi < 2 < -22 7 where = , Z2 = . 



Since —zi < zi, the above inequality then infers that \z\ < Z2. Thus, there exists a 
positive constant C2 = C2(i?, ij) such that < Z2 < Substituting this 

bound into the first term on the left-hand side of (3.7) we also get 



e{n — 1) 



(V''=+i(r))'dr 



ClR 



(3.8) 



Now using the pointwise estimate for linear elliptic equations |16[ Theorem 3.7] 
we have 



max 

[0,-R] 



\^^+\r)\ < ( 



C2 , LfiR) 



(3.9) 



Next, we show that V'r^^ is also uniformly bounded (in k) in [0, R]. To this end, 
integrating (3.2a) over (0, r) after multiplying it by and integrating by parts 

twice in the first term yields 

(n2 - l)[n{n - 1) - 1] 



1 



pn(n— 1) 

1 



gn(n-l)-2^fe+l(^)^^ (3.10) 

(V/(s))""V'+'(s)ds 



er 



1 

n(n — 1) 



for all r G (0, R). Using L'Hopital's rule it is easy to check that the limit as r 0+ 
of each term on the right-hand side of (3.101 is zero. Hence, each term is bounded in 
a neighborhood of r = 0. Moreover, on noting that ip'' > 0, by Schwarz inequality we 
have 



< 



(V''(s))" V'+'(s)rfs 

(rs"("-i)(V''-(s))"''rf5)'( 



(3.11) 



Now in view of ( 3.8 1-( 3.11) we conclude that there exists a positive constant 



C4 = C'4{R,Lf) such that 



max|V'^i(r)|<-g^ 

[0,R] £ 2 



By (3.2a) we get 



-^Lf{r) Vre(0,i?). 



(3.12) 



(3.13) 



Again, using L'Hopital's rule and (2.7 1 it is easy to check that the hniit as r — > 0+ of 
each term on the right-hand side of (3.13) exists, and therefore, each term is bounded 
in a neighborhood of r = 0. Hence, it foUows from (3.91 and (3.121 that there exists 
a positive constant C5 = C5(i?, i/) such that 



(3.14) 



To summarize, we have proved that ||'0'^'''^llc3([o,i?]) ^ R, n,Lf) for j = 0, 1, 2 
and the bounds are independent of k. Clearly, by a simple induction argument we 
conclude that these estimates hold for all fc > 0. 

Step 3: Since ||'0*'llc2([o.i?]) is uniformly bounded in k, then both {ip^}k>o and 
{'0r}fe>o are uniformly equicontinuous. It follows from Arzela-Ascoli compactness 
theorem (cf. [101 page 635]) that there is a subsequence of {'0'^}fc>o (still denoted by 
the same notation) and ij} G C^([0, i?]) such that 

ijj^ — > ip and 4'r — ^ V'r uniformly in every compact set E C (0, i?) as fc 00. 



Testing equation (3.2a) with an arbitrary function x G Cq{{0,R)) yields 



n — 1 



'4'r^\r)xrir) dr + e 

hzT)i^'ir)T''^''-\r)xir)dr= H Lf{r)x{r)dr. 



R 

Setting k 00 and using the Lebesgue Dominated Convergence Theorem, we obtain 



R 



i'r{r)xr{r) dr + e 



R 



n-1 



'ipr{r)x{r) dr 



(3.15) 



+ L „,Mn-i) {W)Tx{r)dr= j Lf{r)x{r)dr. 



Since ip S C^([0, i?]), we are able to integrate by parts in the first term on the left-hand 
side of (3.15), yielding 

eArir) + ' Mr) + ^^„(„_i) (V'(r))" - 2./(r)J x(r) dr = 

for all x e Cq{{0, R)). This then implies that 

-sArir) + '-^^^Mr) + ^^^:^{M)T ' Lf{r) =0 Vr e (0,i?), 



that is, 



-er 



" (^ilV^'^W). + - Lf{r) yr G (0, R). 



Thus, Tp satisfies (2.6a) pointwise in (0, i?). 

Finally, it is clear that ip > in [0, i?], and it follows easily from (3.2b I and 
(3.2c) that ip{0) = i/JriO) — and ipr{R) — £R^~^- So we have demonstrated that 



tp e C^iiO, R]) is a nonnegative classical solution to problem (2.6a)-(2.6c). The proof 
is complete. □ 

Remark 3.1. We note that the a priori estimates derived in the proof are not 
sharp in e. Better estimates will be obtained (and needed) in the next section after 
the positivity of Au'^ is established. 



The above proof together with the uniqueness theorem, Theorem 3.1 and the 
strong maximum principle immediately give the following corollary. 

Corollary 3.3. Suppose r'^'-^f e L^{{Q,R)) and f > a.e. in (0, R), then 
there exists a unique nonnegative classical solution w'^ to problem (2.6). Moreover, 



w" > in {0,R), e C^((0, i?)) if f € C°((0, i?)), and vu" is smooth provided that 
f is smooth. 

Recall that w"^ — r'^~^u^ where and are solutions of (2.4 1 and (|2.6[ ) re- 
spectively. Let be the unique solution to (2.6) as stated in Corollary 3.3 and 
define 



u^r) :=.g(i?) 



1 



— w^s)ds VrG(0,i?). 



(3.16) 



We now show that is a unique monotone increasing classical solution of problem 



(2.4) 



Theorem 3.4. Suppose f e C"((0,i?)) and f > in {0,R), then problem ( [2l | 
has a unique monotone increasing classical solution. Moreover, u'^ is .smooth provided 
that f is smooth. 



Proof. By direct calculations one can easily show that u'^ defined by ( 3.16 1 satisfies 



(2.4). Since > in (0, i?), then u"^ is a monotone increasing function. Hence, the 
existence is shown. 

To show uniqueness, we notice that u*^ is a monotone increasing classical solution 



of problem (2.4) if and only if is a nonnegative classical solution of problem (2.6) 



Hence, the uniqueness of (2.4) follows from the uniqueness of (2.6). □ 



4. Convexity of vanishing moment approximations. The goal of this sec- 
tion is to analyze the convexity of the solution u'^ whose existence is proved in Theorem 



3.4 We shall prove that is strictly convex either in (0,i?) or in (0,R — cqs) for 
some e-independent positive constant Cq. From calculations in Section|2]we know that 
D^u'^ only has two distinct eigenvalues: Ai — u^^. (with multiplicity 1) and A2 = ^uf. 
(with multiplicity n — 1). We have proved that A2 > in (0,i?), so it is necessary to 
show Ai > in (0,i?) or in (0,i? — cqe). In addition, in this section we derive some 
sharp uniform (in e) a priori estimates for the vanishing moment approximations it^ 
which will play an important role not only for establishing the convexity property for 
M^, but also for proving the convergence of in the next section. 
First, we have the following positivity result for Ait^. 

Theorem 4.1. Let u"^ be the unique monotone increasing classical solution of 



problem (2.4) and define := r" m^. Then 

(i) wf. > in (0, R). Consequently, Au"^ > in (0, R), for all e > 0. 

(ii) For any tq G (0, R), there exists an Eq > such that > er"^^ and Au*^ > e 
in (ro,i?) for e G (0,eo). 

Proof. We split the proof into two steps. 

Step 1: Since is monotone increasing and differentiable, we have > in 
[0, i?]. From the derivation of Section [2] we know that r'^^^uf. is the unique 
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nonnegative classical solution of (2.6). Let i/?'^ := w^. By the definition of the Lapla- 
cian A we have 



'At 



(4.1) 



So ip" > in (0, R) infers Au^ > in (0, R). To show ip^ > 0, we differentiate ( |2.6a| ) 
with respect to r to get 



n — 1 



rn—l 



£(n-l)(n-2) (w^)" 



^n{n— 1) 



K)" = r"-V(r) in(0,i?). 



From ( |2.6a| ), we have er" ^(^-s^w^^ 
above two equations yields 



,_\)2+„ (^'^)" ~ r-^/- Combining the 



£(n-l)(n-2) (w^) 



|nn(n — 1) 



Substituting = (p^ into the above equation we get 



e(n-l)(n-2) (w^)' 



j,n{n~l) 

(n- 1)2 



j^(ri-l)2+n / ■ 

(4.2) 



(4.3) 



/ 



> 



since f,Lf,w^ > in (0,i?). This means that vj"^ is a supersolution to a linear 
uniformly elliptic differential operator. By the weak maximum principle we obtain 
(cf. [101 page 329]) 

min^^(r) > mm{0 , ip^O) , ip^R)} = min{0, 0, i?""^e} = 0. 

[0,R] 

Here we have used the fact that ip^{R) = R'^-^Au^R) = R'^-^e. Hence, ip^ > in 
[0,R], so Am^ > in [0,R]. 

It follows from the strong maximum principle (cf. [241 Theorem 4, page 7]) that 
ip"^ cannot attain its nonpositive minimum value at any point in (0, i?). Therefore, 
p'^ > in (0; R), which implies that Au*^ > in (0, R). So assertion (i) holds. 

Step 2: To show (ii), let V'^ '■— w% ~ = r"-^^{Au'^ — e). Using the identities 

— + er"^-^ and w^^^. — + e[n — l){n — 2)r"~^, we rewrite (4.2 ) as 



e{n-l){n-2) [w^) 



e\n—l 



j,7i(n—l) 

(w'^)"-i[(n- l)2w'^ -enr"] 



(4.4) 



in (0,i?). 



Hence, satisfies a linear uniformly elliptic equation. 

Now, on noting that w'^ > by (i), for any tq G (0,i?) (i.e., tq is away from 0), 
it is easy to see that there exists an ei > such that the right-hand side of (4.4) is 
nonnegative in (rp, R) for all e G (0, £i). Hence, is a supersolution in (ro, R) to the 
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uniformly elliptic operator on the right-hand side of (4.4). By the weak maximum 
principle we have (cf. [TUl page 329]) 

min ip%r) > min{0>''(ro), -(/-^(i?)} = min{0, Au''(ro) - e,0}. 

lro,R] 

Again, here we have used the fact that Au^{R} = e. 

Since Au^(ro) > 0, we can choose = minjei, iAu'^(ro)}. We then have 
V^^(ro) = Au^(ro) ~ e > ^Au^ro) > for e G (0,£o). Thus, mm[r„,R]^%r) > 
for e e (0,£o)- Therefore, > er"~^, and consequently, Am^ > e in [ro,i?] for 
£ e (0,eo). 

Finally, an application of the strong maximum principle (cf. [Ml Theorem 4, page 
7]) yields that > er""^. Hence, Au^ > e in (ro,i?) for e € (0,eo). The proof is 
complete. □ 

Remark 4.1. The proof also shows that Eq decreases (resp. increases) as tq 
decreases (resp. increases), andv^ :— Au"^ takes its minimum value e in [rQ,R] at the 
right end of the interval r — R. 

With help of the positivity of Au^, we can derive some better uniform estimates 
(in e) for w'^ and . 

Theorem 4.2. Suppose f e C°((0,i?)) and f > in (0, i?). Let u'^ be the unique 
monotone increasing classical solution to problem (2.4a)-(2.4d|. Define = r"^^u^ 
and v'^ = Au'^ = + ^^^u^. Then there holds the following estimates for sufficiently 
small e > 0: 

(i) \\u'\\c<'{lo,R])+ I \uTdr<Co, 

Jo 

(ii) \\u'\\c^[0,R]) + l|u'1lcn([0,7?]) < Ci, 

(iii) |K||c«([o,i?j) < 

C3 

(iv) ||w^||co([ro,-R]) < n-1 ^0 < ^0 < R, 

er^ 

(v) \\Vr\\c°{[ro..R]) < fa^ip VO < < R, 

fR j-R (-i 

(vi) / K(r)|2dr+ / r2("-i)|t,^(r)pdr < — 



e 

(vii)£ / r-''~'^"°'\v''{r)\'^dr+ I ^{ul{r)Yv'' {r) dr <^ Va < n - 1, 



(viii)£ / r"->^(r)|2dr+ / «(r))"->^(r)p dr < ^ /orn > 3, 



(ix) e\ r^-°'\v'^{r)\^ dr + r^-°'u^^{r)\v' {r)\'^ dr < — for n ^ 2, a < I, 
Jo Jo £ 

where Cj — Cj{R, f,n)>0 for j = 0, 1, 2, • • • ,8 are e -independent positive constants. 
Proof. We divide the proof into five steps 
Step T. Since is monotone increasing, 

maxu^(r) < M^(i?) = g(i?). (4.5) 
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On noting that satisfies equation (2.6a), integrating (2.6a I over (0,i?) and using 
integration by parts on the first term on the left-hand side yields 



-ewl{R) + e{n-l) I -w''[r)dr+- 







f 


f 


dr = 


Jo 


j.n— 1 


Jo 



Lf{r) dr. 



Because w^{R) — eR" and > 0, the above equation and the relation 
imply that 



/' 

Jo 



rU—l 



dr^ K(r)|"dr-<ni?[L/(i?) + £2i?"-2]. 



(4.6) 



It then follows from (|£5|), (|4J|) and ( |3.16| that 

g{R) - nR[Lf{R) + e^R''-'^] " < i/(r) < g{R) Vr e [0, R]. (4.7) 



Hence, is uniformly bounded (in e) in [0, i?], and (i) holds. 



Step 2: Let := Au"" = u% 



;^(r"~i<)^. By \2Aa\ we have 



+ - («)"), = r"- V in (0, R). 



(4.8) 



It was proved in the previous theorem that > e in ( § , i?) for sufficiently small 
e > 0, and it takes its minimum value e a,t r = R. Hence we have v^{R) < 0. We 
note that this is the only place in the proof where we may need to require e to be 
sufficiently small. Integrating (4.8) over (0, i?) yields 



r=R 1 

+ -«)" 
r=0 n 



r=R 



r=0 



LfiR). 



We then have (<(i?))" = nLfiR) + enW^~^v^^{R) < nLf{R), and therefore, 

< (i?) = I < (R) \< (nLfiR))". (4-9) 



Here, we have used boundary condition (2.4c) and the fact that v^{R) < and uj! > 0. 
By the definition of w^{r) :— r"~^uf.{r), we have 



w^R) = \w%R)\< i?"-iK(i?)| < R'-^{nLf{R))" . 
Using the identity v^{r) ~ Au^(r) — u^^i^) + ■^T^"r('')i obtain 

77 — 1 77 — 1 

uUH) = A^i^(i?) - — <.(i?) = s — <(i?). 



Therefore, 



\Urr{R)\ < e + 



n-1 
R 



n — 1 

um\<e+-^inLj{R)y 



(4.10) 



(4.11) 



Step 3: From Theorem 4.1 we have that w^{r) > in (0, i?), and hence, is 
monotone increasing. Consequently, 



maxu;''(r) = max|u;^(r)| < w^iR) < R''^^(nLf(R))" . 
[o,B.] lo,R] I ' 'I V J \ 



(4.12) 
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Evidently, (4.12 ) and the relation 10^(7-) = r" ^u^{r) as well as limr-i-o+ — 

imply that there exists tq > such that 



max zi^(r) = max u^(r) < 



1 



[0,-R] 



[0,-R] ' 



{nLf{R)y. 



(4.13) 



Hence, (ii) holds. In addition, since satisfies the linear elliptic equation (4.2), by 
the pointwise estimate for linear elliptic equations [161 Theorem 3.7], we have 

maxkJ(r)| <£i?"-i + J(i?"-i||/|Uo. +(n-l)||r-iL;|Uoo (4.14) 



+ '^'^^\\r-'K.r\\L^). 

n / 



Since 



, it follows from (4.14) that for any tq > there holds 



max lu'^fr)! = max lAu'^fr) 

[ro,R] [roM] 



(4.15) 



< e 



{-) + ^i3T(i?"-'ll/IU^ + {n-l)\\r-'Lfl 
^tq/ err, \ 



(n- 1)2 



||r-i«)"|U»). 



Thus, (iii) and (iv) are true. 

Integrating ( |4.8[ ) over (0, r) yields 

1 



er"-^< + -«)" = i/ in(0,i?). 



By (4.16) and ( 4.13[ ) we conclude that for any ro > there holds 

Vi I f^y'^"'' V/(-^) 



max k;^(r) = max (Aw^(r))r < - (1 + 



[ro.R] 



lro,R] 



(4.16) 



(4.17) 



Therefore the estimate (v) holds. 



Step 4- Testing (2.6a) with and integrating by parts twice on the first term 

£(« - 1) 



on the left-hand side, we get 



e'R-~'w^{R) + |^^(,)|2dr + i^^K(i?)]^ 



2i? 



^e(n-l), 



2r2 



1 



Combing the above equation and (4.12) we obtain 

rR 



wl.{r)Y dr - 



2r2 



-\w'{r)Ydr + 



nr"("-i) 



|w;^(r)|"+Mr 



Consequently, we have 

|r"-iAM"(r)|2dr + 



e(n-l) 2(„_ 







R 



r"-V;(r)r+idr < i?[£2i?"-2 + Lj:(i?)] {nLf{R)y 
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Hence, (vi) holds. 

Step 5: For any real number a < n ~ 1, testing (4.16) with r^'^v^ and using 
v^{R) — £ we obtain 



-R 



i-i-a , sin- I -a) /"^ „_2- 



r''-^-"\v%r)\' dr 



1 



(4.18) 



+ — «(r))%^(r)dr= -L/(rK(r)dr, 



On noting that > 0, < > 0, and L/(r) = s"^-^ f{r) dr < it follows 

from (|4.18D that 



e(n — 1 — a) 



r"-^-"|?;"(r)|'dr + 



1 /■« 1 



«(r))"w^(r) dr 



(4.19) 



3 nn+3-a|| J-||2 
^ ^jjri-l-a ^ II./ 

~ 2 £n2(n — 1 — a)(n + 3 — a) 



dr Va < n — 1. 



This gives last estimate gives us (vii) 

Ne, 
follows 



Next, recall that v'^ :— Au^ ~ u^^ + ^^w^- Therefore, we can rewrite (4.8) as 



n — 1 



-e(r"-i<)^ + «)"-iw" = r"- V + — — «)" in (0, i?). 
Testing the above equation with r^v'^ for /? > 1 — n and using v^{R) — e, we get 



e'i?"-i+^<(i?) + £ / r"-i+^|z;^(r)|2dr + 
£/3(n + /3 - 2) /'''■^n+/3-3|^e(^)|2^^^ ^ r'' « (r))"->^ (r) ^ dr 



n — 1 



^3^^«+/3-2 



It then follows that 

pR 

-£^R''-^+^v'^{R)+£ r"-^+^\v'^ir)\^dr+^^^ (4.20) 

Jq 2 

_£/3(n+^-2) /■''\.«+;9-3|^e(^)|2^^_^ /''r'3«(r))"->^(r)|2dr 



2 7o 



<-/ r"-i+'3|i;^(r)|2dr+ — / r"^i+'5|/(r)|2 dr 



2£ Jo 



+ {n-l) {u%{T)rv^{r)dT. 
Jo ' 

To continue, we consider the cases n — 2 and n > 2 separately. First, for n > 2, 



it follows from ( |4.19[ ) with a = 1 and ( |4.20| with (3 = that 

r-i? 



e 
2 7o 



R 



r"">^(r)|'dr+ / «(r))"->-'(r)|" dr 



(4.21) 



< 



2e Jo 



r''-^\fir)\'^dr + R" 



2 e{v? - 4) 
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When n ~2,wc note that a = 1 is not aUowed in (4.19). Let a < 1 be fixed in (4.19), 
set /3 = 1 — a in (4.20) we obtain 



(4.22) 







1 



R 



R'Wfl 



e{l — a){5 — a) 



dr 



Hence, (viii) and (ix) hold. The proof is complete. □ 

We now state and prove the following convexity result for the vanishing moment 
approximation u"^ . 

Theorem 4.3. Suppose f e C'^((0,i?)) and there exists a positive constant /q 
such that f > fo on [0, i?]. Let denote the unique monotone increasing classical 
solution to problem ( |2.4a[ )-(2.4d). 

(i) If n — 2,3, then either u"^ is strictly convex in (0, i?) or there exists an e- 
independent positive constant cq such that u'^ is strictly convex in (0, R — cqe). 
(a) If n > 3, then there exists a monotone decreasing sequence {sj}j>o C (0, i?) 
and two corresponding sequences {£j}j>o C (0,1), which is also monotone 
deceasing, and {r*}j>o C (0,i?) satisfying sj \ 0+ as j ^ oo andu^„{sj) > 
and R — r* = 0(e) such that for each j > 0, u'^ is strictly convex in (sj, r*) 
for all e G (0, Ej). 
Proof. We divide the proof into three steps. 

Step 1: Let := r"-i< and := Au^ = u^,, + ^? 



and define rj"^ := r" ^u^^. On noting that r" "^v^ 



be same as before, 
,)^-(n-l)r"-3< = 



Vr 



, the identity (4.8) can be rewritten as 



2e 



If 



'\f + £(3-n)r"-3w^ in(0,i?). 



(4.23) 



So 77^ satisfies a linear uniformly elliptic equation. 

Clearly, ry^(O) = 0. We claim that there exists (at least one) ri € (0, R\ such that 
V'^ifi) > 0. If not, then 77^ < in (0, i?], and u^^. satisfies the same inequality. This 
implies that is monotone decreasing in (0, R]. Since uf.{0) = 0, we have < in 
(0, i?]. But this contradicts with the fact that > in {0,R]. Therefore, the claim 
must be true. 



Due to the factor (3 — n) in the second term on the right-hand side of (4.23), the 
situations for the cases n < 3 and n > 3 are different, and need to be handled slightly 
different. 

Step 2: The case n — 2,3. Since > 0, we have 



eri; 



2e 



> 



in (0,i?). 



(4.24) 



Therefore, 77^ is a supersolution to a linear uniformly elliptic differential operator. 
By the weak maximum principle (cf. |101 page 329]) we have min[o,ri] '7'^('") 5^ 
min{0,?7'^(G),77'^(ri)} = min{0, r;^(ri)} = 0. 

Let r^ = max{ri G (0, i?]; 77^ (ri) > 0}. By the above argument and the definition 
of r* we have 77"^ > in [0, r*], rj^{r^) = if r* 7^ i?, and 77^ < in (r*, i?]. If = R, 
then T^*^ > in [0, i?]. An application of the strong maximum principle to conclude 
that 77^ > in (0,i?). Hence, > in (0,i?). Thus, is strictly convex in (0, i?). 
So the first part of the theorem's assertion is proved. 
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On the other hand, if r^^ < R, we only know that u"^ is strictly convex in {0,r^). 
We now prove that R — = 0(e). By (4.231 and the above setup we have 



in (r*, i?). 



Integrating the above inequality over (r* , r) for r < R and noting that rj^ (r* ) < we 
get —er]^ > ^{r" — r") in (r*,i?). Integrating again over (r*,i?) and using the fact 

that r]'^ (r* ) = and the algebraic inequality ^ „ r" > 

at 



71+1 



i?", we arrive 



n(n + 1) 



It follows from (|4.11|) that 
R 



en(n + l)\urr{R)\ en{n + l) 
< r^^^ < 



Rfo - R^fo 

Thus, R — r^, ~ 0{e), and is strictly convex in (0, R 



Cos. 



CqE 



Step 3: The case n > i: First, By the argument used in Step i, it is easy to show 
that can not be strictly negative in the whole of any neighborhood of r = 0. Thus, 
there exists a monotone decreasing sequence {sj"}j>o C (0,i?) such that Sj \ 0^ as 
j — )■ c» and tfisj) > 0. 

Second, we note that 



r(3 — n)r 



e(3- 



„«-3 



1 



e(3-n) ( \ + (n-l)r 



Using this identity in (4.231, we have 



(n-l)e {ulY 



T] = r 



^[r3/ + £(3-n)(n- !)<] in(0,i?). 



By (ii) of Theorem 4.2 we know that is uniformly bounded in [0, R]. Then for each 
Sj there exists an Sj > (without loss of the generality, choose Sj < ej-i) such that 
for £ G (0,£j), there holds [r^f + £(3 — n){n — l)uf.] > in {sj,R). Hence, rj^ is a 
supersolution to a linear uniformly elliptic operator on (sj, R) for £ < £j. 

Third, for each fixed j > 1, let r* = max{r € {sj, R\; 77^ (r) > 0}. Trivially, by the 
construction, r* > Sj-i > Sj. By the weak maximum principle (cf. ^IQi page 329]) 
we have 

min ?7^(r) > minjO, 77'^(sj), ?f (r*)} > 0. 
Finally, repeating the argument of Step 2:, we conclude that is either strictly 



convex in {sj,R) or in {sj,r*) with R 
complete. □ 



0{e) for £ £ (0,£j). The proof is now 



5. Convergence of vanishing moment approximations. The goal of this 
section is to show that the solution of problem (2.4) converges to the convex 



solution u of problem (2.1 ) We present two different proofs for the convergence. The 



first proof is based on the variational formulations of both problems. The second 
proof, which can be extended to more general non-radially symmetric case, is done 
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in the viscosity solution setting [9|. Both proofs mainly rely on two key ingredients. 
The first is the solution estimates obtained in Theorem |4.2[ and the second is the 
uniqueness of solutions to problem ( |2.1[ ). 

Theorem 5.1. Suppose f e C"^((0, R)) and there exists a positive constant fo 
such that f > fo in [0, i?]. Let u denote the convex (classical) solution to problem 
(2.1) and he the monotone increasing classical solution to problem (2.4 1. Then 

(i) — limg_>Q+ exists pointwise and converges to uniformly in every 
compact subset of (0, i?) as e . Moreover, vP is strictly convex in every 
compact subset. Hence, it is convex in [0,i?]. 

(ii) u^j, converges to weakly * in L°°{{0, R)) as e ^ 0+. 

(iii) = u. 



Proof. It follows from (ii) of Theorem 4.2 that ||u^||(71([o,_r]) is uniformly bounded 
in e, then {it'^}£>o is uniformly equicontinuous. By Arzela-Ascoli compactness theo- 
rem (cf. [ini page 635]) we conclude that there exists a subsequence of {w^}e>o (still 
denoted by the same notation) and G C^([0, i?]) such that 



u 

ul 



uniformly in every compact set E C (0, R) as £ — ^ O"*", 
weakly * in L°°((0, i?)) as e -^■ 0+, 



and u'^{R) = g{R) implies that 



u°{R) = g{R). 



In addition, by Theorem 4.3 we have that for every compact subset E C (0, R) 
there exists eg > such that E C {0,R — cqe) and is strictly convex in (0, R — cqe) 
for e < Eq. It follows from a well-known property of convex functions (cf. [5Dj) that 
must be strictly convex in E and S Cl,'}{{0, R)). 

Testing equation (4.8 1 with an arbitrary function x G Cq{{0,R)) yields 



r" '-v^{r)xr{r) dr - 



{ul{r)Yxr{r)dr= / r'^-^ f{r)x{r) dr, (5.1) 



where as before — Au^ — u^^-t-^^y^w^. By Schwartz inequality and (vi) of Theorem 
4.21 we have 



r" 'vl{r)xr{r)dr = -e / r" 'v^{r) Xrr{r) + 



1 



< e 



< \/^5 



r^(^-^)\v^{r)\^dr 
\Ax{r)\^dr 



r 

R. 



-Xr{r) 



|Ax(r 
as e - 



dr 



l^dr 



0^ 



Setting e 0+ in (5.1 1 and using the Lebesgue Dominated Convergence Theorem 
yields 



{ul{r)Yxr{r)dr= r"-' f{r)x{r) dr Vx G Ci((0, i?)). 



(5.2) 



It also follows from a standard test function argument that u"(0) = 0. This means 
that e C"'^([0, R]) n Ci'<-^((0, R)) is a convex weak solution to problem (2.1 1. By the 
uniqueness of convex solutions of problem (2.1 ), there must hold u'^ = u. 
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Finally, since we have proved that every convergent subsequence of {u'^}e>o con- 
verges to the unique convex classical solution u of problem (2.1), the whole sequence 
{u^}e>o must converge to u. The proof is complete. □ 

Next, we state and prove a different version of Theorem |5.1[ The difference is 
that we now only assume problem (2.1 ) has a unique strictly convex viscosity solution 
and so the proof must be adapted to the viscosity solution framework. 

Theorem 5.2. Suppose f G C'^{{0, R)) and there exists a positive constant /o 
such that f > fo on [0, R] . Let u denote the strictly convex viscosity solution to 
problem (2.11 and he the monotone increasing classical solution to problem (2.4). 
Then the statements of Theorem |5.1| still hold. 



Proof. Except the step o f pro ving the variational formulation (5.2), all other 
parts of the proof of Theorem 5J_ are still valid. So we only need to show that vP 
is a viscosity solution of problem (2.1), which is verified below by the definition of 
viscosity solutions. 

Let e C'^([0, i?]) be strictly convex. Suppose that ~ (f> has a local maximum 
at a point ro € (0,i?), that is, there exists a (small) number 5q > Q such that 
(ro - i5o, ro + <5o) CC (0, R) and 

u\r) - (j){r) < u°(ro) - 0(ro) Vr € (ro - 5o,ro + So). 

Since (which still denotes a subsequence) converges to uniformly in [ro — 
So, ro + So], then for sufficiently small e > 0, there exists r^ G (0, R) such that r^ — > ro 
as e — i' 0+ and — (f) has a local maximum at r^ (see [TUl Chapter 10] for a proof of 
the claim). By elementary calculus, we have u%{r^) = (t>r{re) and u'^„{r^) < (f>rr{re)- 
Because both and <f) are strictly convex, there exists a (small) constant po > such 
that for sufficiently small £ > 

ulrir) < (t)rr{r) Vr e (ro - p, ro + p), p< po, 

which together with an application of Taylor's formula implies that 

<(r) = (j)r{r) + 0{\r - r^\) Vr E (ro - p,ro + p), p < po. 



Let X e C'o((''o-/0, ro + p)) with x > and x(ro) > 0. Testing (4.8) with x yields 



1 

2np 



ro+P 



> 



ro-p 
1 

2p 
1 

2np 
1 



2p 



ro+P 



ra~p 



.(r))"-Vrr(r)x(r)dr 



ro+P 

[«(r))"-i + 0(|r - r,|"-i)]<,(r)x(r) dr 

ro-p 

ro+P r -, 

(«(r))")^ + 0(|r - r,r-i)<,(r) x{r) dr 



> 



2p 



ro-p 
ro+P 



r''-' [/(r)x(r) + e<(r)xr(r)] dr - C^p^- 



ro+P 



(5.3) 



K['r)Xr{r) dr 



for some positive p- independent constant Cg. Here, we have used the fact that > 
0, X ^ in [ro — p, ro + p] to get the last inequality. 

18 



From (vi) of Theorem 4.2 we have 



ro+p 



r^-^vl{r)xr{r)dr 



< e 



ro+p p n - 1 1 

r«-i„^(r) xrr{r) + Xr{r) dr (5.4) 

T 

ro-p 
R 



< 



'^\v%r)\'dry(^J^ ' \Ax{r)\-' dry 
^,[r \Ax{r)? dry . 



Setting e — > 0+ in (5.3) and using (5.4) we get 

Cro+p 



2p 



.(r))"-Vr.(r)x(r) dr 



(5.5) 



> 



2p 



ra+p 



r''-'fir)xir)dr-Cgp'' 



TQ-P 



ra+p 



Ur{r)Xr{r) dr. 



Here, we have used the fact that uf. converges to m° weakly * in L°°{{0,R)) to pass 
to the hmit in the last term on the right-hand side. 

Finally, letting p — > 0+ in ( |5.5[ ) and using the Lebesgue-Besicovitch Differentiation 
Theorem (cf. [TO]) we have ((?!)r(?'o))"^ Vrr(?'o)x(»^o) > r^'^ f{ro)x{ro)- Hence, 



0rr(fo) > /(?'o), 



so is a viscosity subsolution to equation (2.1a). 

Similarly, we can show that if — assumes a local minimum at ro G (0, R) for 
a strictly convex function (p G C!q{{0, R)), there holds 

-(/)^(ro)n"-l 



ro 



Therefore, m° is also a viscosity supersolution to equation (2.1a). Thus, it is a viscosity 
solution. The proof is complete. □ 

6. Rates of convergence. In this section, we derive rates of convergence for 
u'^ in various norms. Here we consider two cases, namely, the n-dimensional radially 
symmetric case and the general n-dimensional case, under different assumptions. In 
both cases, the linearization of the Monge-Ampere operator are explicitly exploited, 
and it plays a key role in our proofs. 

Theorem 6.1. Let u denote the strictly convex classical solution to problem 
(2.1) and he the monotone increasing classical solution to problem (2.4). Then 
there holds the following estimates: 

(6.1) 
(6.2) 

where Cj = Cj(||r' 
and 



e'ir)\urir)~uUr)\'dr) < do, 



\Au(r) - Au^r)f dr) < Cn 







' ^Awrllis) for j = 10, 11 are two positive e-independent constants, 
^'''•^^ - Y.^Ur{r)Y{ul{r)r-^-= > in(0,i?]. (6.3) 
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Proof. Let v :— Au — 



^-Ur, := Au"^ = u% — 



„ ^^u^, and :— u — 



On noting that (2.4a) can be written into (4.81, multiplying (2. la I by r" and 
subtracting the resulting equation from (4.81 yields the following error equation: 



1 



(r-ix;^)^ + - [{urT - (O'l - in (0, R). 



Testing (6.4| with e"^ using boundary condition e^(0) = e^{R) = we get 



1 



^(r)|e^(r)prfr = 0, 



(6.4) 



(6.5) 



where 9'^ is defined by (6.3 1. 



Integrating by parts on the first term of (6.5 1 and rearranging terms we obtain 
rR 



fR . fR 

/ r"-i|Ae"(r)p(ir+- / ( 
J a ^ Jo 



e%r)\eUr)\' dr 



(6.6) 



eR''-'Ae'{R)e^{R) ~ e / r"-^w^(r)e^(r) dr. 



R 



We now bound the two terms on the right-hand side as follows. First, for the 
second term, a simple application of the Schwarz and Young's inequalities gives 



r''-'vr{r)e';.{r)dr < 



1 

4n Jo 



(r)|e^(r)|^dr + £^n 



R j,2{n~l) 

dHr) 



\vj.{r)\'^ dr. 



Second, to bound the first term on the right-hand side of (6.6), we use the boundary 
condition v'^{R) = e to get 



and 



|Ae"(i?)| = \v{R) ~ v'{R)\ = \v{R) - e\ < \v{R)\ + 1=:M, 

pR. pR 

\R"-^e'^{R)f ^ {{r'^-^el{r)f)^dr ^2 r2("-i)e^(r) Ae^(r) dr 
^0 ^ Jo 

.l.j-R ,1 
r"-i|Ae^(r)|2 dry U r3("-i)|e^(r)p drj ' . 



< 2 



Hence by Young's inequality, we obtain 



ati. 1 pR. 

r"-i|Ae^(r)|2 dr) ' r3("-i)|e^(r)|2 dr 



< 



< 



r''-^\Ae'{r)\'^dr + 2sMU / r^^"-!) |e^(r)|2 dr 



r"-^|Ae'^(r)|^dr + 



4n 



(6.7) 



0%r)\e'^{r)\'dr + e^nM'C 



for some e-independent constant C = C{f, R, n) > 0. 
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Combining (6.6)-(6.7) yields 



r"-i|Ae^(r)|2dr + 



'irMir)\'dr 



(6.8) 



< 2e^n 



R „2(n-l) 



-\vr{r)\'^ dr + e'inM'^C. 



Thus, (6.1) and (6.2) follow from the fact that \\r 



L°o < OO. □ 



Corollary 6.2. Inequality (6.1) implies that there exists an e-independent con- 
stant C > such that 



R 



r"'-^\ur{r)-ul{r)\'Ur]^ <e^CC_ 



-•lO- 



(6.9) 



Since the proof is simple, we omit it. 



Theorem 6.3. Under the assumptions of Theorem 6.1, there also holds the 
following estimate: 



. i-R ,1 

(/ r"-^\u{r)-u%r)\^dry <eCi 



(6.10) 



for some positive e-independent constant C12 — Ci2(i?, n, u, Cn). 



Proof. Let 0^ be defined by (6.3), and e^, v and be same as in Theorem 6.1 



Consider the following auxiliary problem: 



nr'^-^e" in(0,i?). 



(6.11a) 
(6.11b) 



We note that the left-hand side of (6.11a) is the linearization of (2.1al at 9' 



Since 6*^ > in (0,i?], then (6.11a) is a linear elliptic equation. Using the fact 
that ci > r"~^(6'^)~^ > cq > in [0,i?] for some e-independent positive constants 
Co and Ci, it is easy to check that problem (6.11) has a unique classical solution <f>. 
Moreover, 

R pR fR. 

r''-'^\(j)rr{r)\'^ dr + r"-i|</)^(r)p dr < C / r"-i|e^(r)|2 dr (6.12) 
Jo Jo 

for some e-independent constant C = C(/, i?, n, cq, ci) > 0. 

Testing (6.11a) by e^, using the facts that 0^(0) — i'iR) — 0, e^(i?) — and 
v^{R) = e as well as error equation \6A\ we get 



R 1 j-R i-R 

r""'^\e'{r)\'^dr^ / ^ (r)0^(r)e^(r) dr = e / r''-\'^{r)<Pr{r) dr (6.13) 

"Jo Jo 

fR 

= eR'''^v'{R)(f>r{R)-e / r''-^v'{r)A(j){r) dr 
Jo 

I-R 

= e^i?"- Vr(i?) +s r"-i [v(r) - v%r)]A(j){r) dr 
Jo 

.R 

-e r''-^v{r)A(j){r)dr, 
Jo 
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where we have used the short-hand notation = r 



n-l r 



For each term on the right-hand side of (6.13) we have the fohowing estimates: 



r"-i[u(r) - w^(r)]A0(r)dr 



dr 



< £ 



1 , nR 



''-^\v{r)~v'{r)fdry U r"-i|A(/.(r)|2 dr) 







1 , pR 



r''-\{r)A4>{r)dr<e(^J r"-Xr)|2 dr) ' r"-i|A(/.(r)|2 dr) 



Substituting the above estimates into (6.13) and using (6.12) we get 

cR 

r"-i|e''(r)pdr 



(6.14) 







(/ r-"-i|0,,(r)|2dr)' + ^(/ r^-'\Mr)\^ dr) 



r'^-^\A(j){r)Y dr 

R 



< Ae{eR^ +eV^R-^ +e^Cu+Cu)c(^J^ r"-i|e''(r)|2 dr 

for some e- independent constant C„ — C{u) > 0. 

Hence, by ( 6.14| we conclude that ( |6.10| ) holds with C12 = 4(ei?t + £y/nR~'^ + 

£^C\\ -I- Cu)C. The proof is complete. □ 

Since the proofs of Theorem |6.1| and |6.3| only rely on the ellipticity of the lin- 
earization of the Monge- Ampere operator, hence, the results of both theorems can be 
easily extended to the general Monge- Ampere problem ( 1.3 ) and its vanishing moment 
approximation (1.2). 

Theorem 6.4. Let u denote the strictly convex viscosity solution to problem ( 1.3 ) 
and be a classical solution to problem (1.2) Assume u £ W'^'°^{Q.) r\H^(Q?) and 
is either convex or "almost convex" in (which means that is convex in Q minus 
an e -neighborhood of dQ; see Theorem 4.3 for a precise description). Then there holds 
the following estimates: 



Vu-Vu^pdxj' <eiCi3, 
1 

Au - A-U^P dx) ' < £2 Cl4, 

u - u^p dx) ^ < eCi5, 



(6.15) 
(6.16) 
(6.17) 
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where Cj ~ Cj{\\V Au\\l2) for j = 13, 14, 15 are positive e- independent constants. 
Proof. Since the proof follows the exact same lines as those for Theorem |6.1[ we 



just briefly highlight the main steps. First, the error equation (6.4 1 is replaced by 
eAt;^ + del{D'^u) - Aet{D'^u'') = 

where v 



in 57, 



(6.18) 



e / lAe' 
In 



Au^ . Next, equation (|6.3|) becomes 
'^^dx + 



'Ve" ■ Ve" dx 



n 



Ae"— dS^- 



e I VvVe^dx, (6.19) 
In 



where 



= <^>^ := cof(tD'^u + (1 - t)D^u''^ 



for some t G [0,1], 



(6.20) 



now stands for the cofactor matrix of tD^u + (1 — t)D^u^ . Since u is assumed to be 
strictly convex, and it^ is "almost convex" , there exists a positive constant Of) such 
that rVe^ • Ve^ > 6'o|Ve=p. 



It remains to derive a boundary estimate that is analogous to (6.7 1. To the end 



by the boundary condition v'^\qq = e and the trace inequality, we have 



Ae'—ds < e(e\dn\ 
dn dv 



\^Al^on)) 



de' 



dv 



L^(dn) 



(6.21) 



< eM||Ve^ 



L2(r!) 
< |l|Ae1li2(o)-l- 



lAe 



M3£||Ve= 



\h{n) 



< 



|Ae^"2 



The desired estimates (6.15) and (6.16) follow from combining (6.19) and (6.21) 



Finally, (6.17) can be derived by using the same duality argument as that used 



in the proof of Theorem |6.3[ We leave the details to the interested reader. □ 

7. Numerical Experiments. In this section, we numerically verify the theo- 



retical results given in the previous sections. To this end, we solve (2.4a)-(2.4d) in 
the domain = (0, 1). We use the Hermite cubic finite element to construct our finite 
element space (cf. [3]), and we use the following data: 



f = il + r') 



5(1) 



62 . 



It can be readily checked that the exact solution is w = e' 



V2 



We plot the computed solution and corresponding error in Figure 7.1 with pa- 
rameters n = 4, e = 10~^,h — 4.0 x 10~^. We also plot the computed Laplacian, 
Au^ := u^j, + ^M^, as well. As shown by the pictures, the vanishing moment methodol- 
ogy accurately captures the convex solution in higher dimensions. Also, as expected. 



the Laplacian of u'^ is strictly positive (cf. Theorem 4.1). 

Next, we plot both Ur and Urr in two and four dimensions in Figures 7.2-7.3 with 
£-values, 10^^, 10^"^, 10^^. Recall that the Hess ian matrix of u'^ only has two distinct 
eigenvalues u^^ and ^u^. As seen in Figure 
both dimensions n = 2 and n 



7.2 



Finally, Figure 7.3 shows that u: 



is positive for all e-values and for 
4. This result is in accordance with Corollary |3.3[ 
is strictly positive except for a small e-neighborhood 



of the boundary, which agrees with the theoretical results established in Theorem 4.3 
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Computed Solution: n=4, e=lE-3 



Laplacian; n-4. e-IE-3 



0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Fig. 7.1. Computed solution of ( |2.4| (left) and and computed Laplacian (right) with n = 4,£ : 
10-1. h = 4 X 10-^. 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Fig. 7.2. Computed Ur of \2A\ for n = 2 (left), and n = 4 (right) with e = 10~^ (black), 
e = 10-3 (blue), and e = IQ-^ (redf(h = 4 X 10~''^). 



8. Concluding Remarks. We like to comment on an interesting property of 
the vanishing moment method, namely, the ability of the vanishing moment method 



to approximate the concave solution of the Monge- Ampere problem (1.3). This can 
be achieved simply by letting e 0~ in (1.2 1. This property can be easily proved as 
follows in the radially symmetric case. 



Before giving the proof, we note that for a given / > in il, equation ( 1.3a) does 
not a have concave solution in odd dimensions (i.e., n is odd) because det{D^) = f 
does not hold for any concave function u as all n eigenvalues of Hessian D^u of a 
concave function u must be nonpositive. On the other hand, in even dimensions (i.e.. 



n is even), it is trivial to check that if u is a convex solution of problem (1.31 with 
0, then —u, which is a concave function, must also be a solution of problem (1.3). 



9 



Next, by the same token, it is easy to prove that if is a convex or "almost 

, which is concave or "almost concave" 

in the sense that 



convex" solution to problem (1.2), then - 

-u^ is concave in 57 minus an 0(e)-neighborhood of the boundary 



dfl of f2), must also be a solution of (1.2) 



Finally, let n be a positive even integer, it is easy to see that changing u"^ to — 
in (1.2) is equivalent to changing e to —e in (1.2). For e < 0, let 5 := — e. After 
replacing e by —5 and by := — in (2.4 1, we see that satisfies the same 
set of equations (2.4) with 6{> 0) in place of e. Hence, by the analysis of Sections 
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6| we know that there exists a monotone increasing solution to problem (|2.4[) 
with e being replaced by S, which satisfies all the properties proved in Sections ^||6j 
Translating all these to = —u^, we conclude that problem (2.4 1 for e < has a 
monotone decreasing solution which is either concave or "almost concave" in (0, R) 
and converges to the unique concave solution of problem ( 1.3 ) as e 0^ . In addition, 
u"^ satisfies the error estimates stated in Theorems 16.11 and 16.31 
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